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Abstract. Program SHELLSPEC is designed to calculate light-curves, spectra and images of 
C^i interacting binaries and extrasolar planets immersed in a moving circumstellar environment 
which is optically thin. It solves simple radiative transfer along the line of sight in moving 
^jpj media. The assumptions include LTE and optional known state quantities and velocity fields in 
r-{ I 3D. Optional (non)transparent objects such as a spot, disc, stream, jet, ufo, shell or stars may 
be defined (embedded) in 3D and their composite synthetic spectrum calculated. Roche model 
can be used as a boundary condition for the radiative transfer. Recently a new model of the 
lO ] reflection effect, dust and Mie scattering were incorporated into the code. 

I e Aurigae is one of the most mysterious objects on the sky. Prior modeling of its light- 

curve assumed dark, inclined, disk of dust with the central hole to explain the light-curve 

with a sharp mid-eclipse brightening. Our model consists of two geometrically thick flared disks. 
Internal optically thick disk and external optically thin disk which absorbs and scatters radiation. 
Shallow mid-eclipse brightening may result from eclipses by nearly edge-on flared (dusty or 
, gaseous) disks. Mid-eclipse brightening may also be due to strong forward scattering and optical 
Oh' properties of the dust which can have an important effect on the light-curves. 
^ I There are many similarities between interacting binary stars and transiting extrasolar planets. 

5^ Reflection effect which is briefly reviewed is one of them. The exact Roche shape and temperature 
distributions over the surface of all currently known transiting extrasolar planets have been 
■ determined. In some cases (HAT-P-32b, WASP-12b, WASP-19b) departures from the spherical 
' ' shape can reach 7-15%. 

T-H Keywords, binaries: eclipsing, planets and satellites: general 
> 

5^ ' 1. Introduction 

There are sophisticated computer codes for calculating and inverting light curves or 
r-^ spectra of binary stars with various shapes or geometry including the Roche model (Lucy 
^ ; 1968; Wilson & Devinney 1971; Wood 1971; Mochnacki & Doughty 1972; Rucinski 1973; 

I HiU 1979; Popper & Etzel 1981; Djurasevic 1992; Drechsel et al. 1994; Hadrava 1997; 
^ Bradstreet & Steelman 2002; Pribulla 2004, Pavlovski et al. 2006, Tamuz et al. 2006). 
The Wilson & Devinney code is most often used and is continuously being improved or 
modified (Kallrath et al.l998; Prsa & Zwitter 2005). The main focus of these codes is 
^ to deal with the stars, determine their properties and their orbit. However, it is often 
" " ' the case that stellar objects are embedded in some moving optically thin environment 
and/or are accompanied by discs, streams, jets or shells which give rise to various emission 
spectra. We address these objects with Shellspec. 

The Shellspec code is described in more detail in Budaj & Richards (2004). There has 
been a lot of progress since that manual was published and an updated version of the 
manual, with examples of input, output and test runs is available within each new release. 
A convenient overview may also be found in Budaj & Richards (2010). Major changes 
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since that time include a new model of the reflection effect applicable to tidally distorted 
and strongly irradiated cold objects (Budaj 2011a), and dust and angle dependent Mie 
scattering (Budaj 2011b). It includes extinction due to the absorption and scattering as 
well as the thermal and scattering emission. Original Shellspec code is written in For- 
tran77 and does not solve the inverse problem of finding the best fit parameters. There 
are versions of this code in Fortran90 which even solve some particular restricted inverse 



problems by Tkachenko et al. (2010) and Sejnova et al. (2011). Chadima et al. (2011a^ 



modeled Ha emission V/R variations caused by discontinuous mass transfer in interact- 
ing binaries. Miller et al. p007|) modeled UV and optical spectra of TT Hya, Algol type 
binary. Sections below describe applications of this code to e Aur and extrasolar planets. 
Output of the codes SYNSPEC and COOL-TLUSTY (Hubeny 1988; Hubeny, Lanz & 
Jeffery 1994; Hubeny, Burrows & Sudarsky 2003) is used as a default input of spectra of 
nontransparent objects, which serves as a boundary condition for the radiative transfer 
in the interstellar medium. 



2. e Aurigae 

e Aur is an eclipsing binary with the longest known orbital period, 27.1 yr. A very rare 
eclipse that lasted for two years is over. However, the object, its origin and, in particular, 
the secondary component of this binary star remain mysterious. Huang (I1965P proposed 
that the secondary is a dark disk seen edge-on. Wilson (|197ip and Carroll et al. (|199ip 
argued that the observed sharp mid-eclipse brightening (MEB) can only be explained 
by a tilted disk with a central opening. Ferluga ((1990,) suggested that the disk is a system 
of rings. 
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Figure 1. Eclipse of e Aur by a dark, geometrically thick, flared disk of dust. The disk consists 
of two parts; (1) The flared optically thick part that causes most of the eclipse, and (2) a flared 
optically thin part that causes additional absorption, scattering and mid-eclipse brightening. 
Model (A): Disk has only one part (1). Model (B): disk has both parts part (1) and part (2). 
Mid-eclipse brightening arises mainly because the edges of the flared disk are more effective in 
the attenuation of the stellar light than the central parts of the disk. Thin dotted line is a best-fit 
quadratic function to the eclipse bed. Crosses - observations from AAVSO (Henden l201ip . 

There has been a wealth of studies during the current eclipse. Orbital solutions were 

f By a mid-eclipse brightening in e Aur some authors understand only a relatively sharp local 
maximum that appeared in a few eclipses near the middle of the eclipse. We propose a slightly 
more general definition of the MEB or a new term (mid-eclipse excess). Our MEB or mid-eclipse 
excess is a convex feature near the middle of an eclipse bed. A common eclipse has a concave 
eclipse bed. 
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Figure 2. 2D image of e Aur during eclipse. Black region is a dark, geometrically and optically 
thick, edge-on flared disk of dust. This disk causes the most of the eclipse. Colored regions 
correspond to the optically thin flared disk which scatters and absorbs the light from the F-star. 
Notice that this optically thin disk may produce artificial spots on the surface of the F star. 

recently revisited by Stefanik et al. (|2010p and Chadima et al. pOlOp . Kloppenborg et al. 
(j20101 1201ip confirmed the dark disk with interferometric observations but they did not 
confirm the hole in the disk. The spectral energy distribution was studied by Hoard et 
al. (PITUI) . These authors favor the post AGB+B5V model. Wolk et al. (POTU)) analyzed 
X-ray observations. Sadakane et al. (|2010! ) carried out the abundance analysis. Recently, 
Chadima et al. p010[ I2011bp questioned the presence of sharp mid-eclipse brightening 
and suggested that the photometric variability seen during eclipse is intrinsic to the F- 
star. Extremely precise observation of the light-curve of e Aur were obtained with the 
Solar Mass Ejection Imager by Clover et al. (2011). These data clearly show a shallow 
mid eclipse brightening. 

Observations that were used for comparison with our calculations were taken from the 
AAVSO database (Henden l201ip . They were obtained by many observers who contributed 
to the database during the current eclipse of e Aur. Only the observations in the V-filter 
were considered here. These observations also indicate the presence of a shallow mid- 
eclipse brightening but it is not as sharp and pronounced as some might have anticipated. 
Ingress is not as steep as egress, which indicates that the disk is not perfectly symmetric, 
but suffers from some small disturbance. Its leading part might be more extended or disk 
slightly inclined out of the orbital plane (warped?). 

To explain these observations we suggested an alternative model of e Aur. Our model 
of e Aur consists of two geometrically thick flared disks: an internal optically thick disk 
and an external optically thin disk, which absorbs and scatters radiation. Disks are in 
the orbital plane and are almost edge-on. The model is based on optical properties of 
dust grains. It takes into account the extinction due to the Mie scattering and absorption 
as well as thermal and scattering emission. We argue that there is no need for a highly 
inclined disk with a hole to explain the current eclipse of e Aur even if there is a possible 
shallow mid-eclipse brightening (see Fig[T]). For mode details kindly see Budaj (2011b). 
Figl2] displays a 2D image of e Aur as calculated in the V band. 

3. Shapes of transiting extrasolar planets 

Transiting extrasolar planets are very close to their parent stars. Most of them have 
circular orbits or very small eccentricity. This indicates that their rotation is synchronous 
and classical Roche model can be applied to them. Budaj (2011a) calculated the shape 
of all transiting exoplanets known at that time. By the shape we mean the relative 
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Table 1. Shapes of the transiting exoplanets. Columns are: a -semi- major axis in [AU], Rsub 
-planet radius at the sub-stellar point, Rback -planet radius at the anti-stellar point, RpoU 
-planet radius at the rotation pole, Rside -planet radius at the side point, (assumed equal to the 
planet radius determined from the transit), Reff -effective radius of the planet, rr = Rsub/Rpou 
-departure from the sphere, fi = Rsut/Lix -fill-in parameter of the Roche lobe. Radii are in 
units of Jupiter radius. See the text for a more detailed information. 

a Rsub Rback Rpole Rside RefF rr fi name 



0.01655 1.55474 1.54069 

0.02293 1.90438 1.89283 

0.03440 2.14496 2.14005 

0.02540 1.54395 1.54167 

0.02312 1.40758 1.40572 

0.05150 2.04439 2.04295 

0.02250 1.23086 1.22972 

0.03444 1.70957 1.70825 



1.35022 1.38600 1.42045 

1.69728 1.73600 1.77110 

2.00804 2.03700 2.06019 

1.47462 1.49000 1.50166 

1.35013 1.36300 1.37264 

1.97515 1.99100 2.00266 

1.19084 1.20000 1.20670 

1.65814 1.67000 1.67862 



1.151 0.63 WASP-19b 

1.122 0.58 WASP-12 b 

1.068 0.48 HAT-P-32 b 

1.047 0.42 CoRoT-1 b 

1.042 0.41 WASP-4 b 

1.035 0.38 WASP-17b 

1.033 0.37 OGLE-TR-56 b 

1.031 0.36 WASP-48b 



proportions of the object. Shape does not change much if the absolute dimensions change 
within typical measurement errors. Consequently, the new absolute dimensions can be 
easily linearly rescaled e.g. if the measured cross-section of the planet changes. I use 
Rsub /Rpole as a measure of the departure from the spherical shape, where Rsub is the 
radius at the sub-stellar point and Rpoie is the radius at the rotational pole. I updated 
these calculations and recalculated the shape of all 138 transiting exoplanets known as 
of end of June, 2011. Most of the transiting exoplanets have small departures from the 
sphere of the order of 1%. However, a fraction of them have departures more than 3% and 
these are Usted in the Table 1. HAT-P-32 b, WASP-12 b, and WASP-19 b are recorders 
with departures of about 7. 12, and 15% respectively. This means that we are observing 
only a cross section of these planets during the transit. This cross-section is also not 
spherical, it is characterized by the Rside /Rpole parameter but this is not as high as the 
Rsub/ Rpole- Table 1 lists also the effective radius of the planet (radius of the sphere with 
the same volume) which one could use e.g. for comparison with the theoretical radius 
calculations. Table 1 illustrates that spherical shape of close-in exoplanets cannot be 
taken for granted. Moreover, the Table 1 demonstrates that this Roche shape cannot be 
approximated by the rotational ellipsoid. This assumption used in some papers is not 
justified for most of the transiting extrasolar planets. 



4. Reflection eflfect 

Reflection effect operates in very different environments and there are different ap- 
proaches to model this effect. The table below presents an overview of how the models 
work in the field of interacting binaries and exoplanets. One can see that the models 
of the reflection effect and definition of the fundamental quantities are very different in 
these fields. The question is: Is there any transition between the two approaches? What is 
the amount of the heat redistribution? What fraction of light gets reflected off the surface 
and does not produce heat? Which approach (model) should one choose for tidaly dis- 
torted and strongly irradiated cold object? To address these questions and these objects 
Budaj (2011a) proposed a sort of hybrid model of the reflection effect whose properties 
are also listed in the table. 
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Table 2. Overview of the reflection efiiect. Scattering refers to the refletion of Ught off the 
surface of one of the objects. Aboi is bolometric albedo (Rucinski 1969), Ab is Bond albedo, Av 
is monochromatic albedo. MR stands for muliple reflection between surfaces of the two objects. 



Property Interacting Binaries Exoplanets Hybrid model 

Shape Roche (limb+grav.dark.) sphere, rot. ellipsoid Roche (limb+grav.dark.) 

Scattering No Yes Yes 

Albedo Aboi Ab,A^ Ab,A^ 

(A) Meaning absorbed=>heating reflected=>no heating reflected=>no heating 

(1-A) Meaning penetrates into star heating+heat redistrib. heating+heat redistrib. 

MR Yes No No 
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Discussion 

Wilson: Did you look into the possible importance of multiple scattering? The scat- 
tering geometry affects the angular distribution of scattered light and, to a small extent, 
even whether a scattered photon escapes the atmosphere. The effect is small but might 
not be negligible. 

BuDAJ: There is a scattering (reflection) process in two environments. One is the scat- 
tering in the atmosphere of stars or nontransparent objects. I do not solve the radiative 
transfer in the atmospheres of these objects. I take the flux from the model atmosphere 
calculation and use is to calculate the boundary condition for the radiative transfer in 
the moving interstellar matter. Models of the atmosphere take the multiple scattering 
into account. The other is the scattering process in the interstellar medium. In this case I 
assume that the medium is optically thin and that it is irradiated by one or two sources. 
I take into account only the first scattering event. In the optically thin medium the prob- 
ability of the second and higher scattering events rapidly decreases. To take them into 
account one would have to solve the 3D radiative transfer in the 3D moving medium. 

Burrows: Have you thought of using your code to calculate light curves for WASP-12b 
at the warm Spitzer bands (3.6 m and 4.5 m), instead of at 8 m? There should be some 
data soon from Spitzer for WASP- 12b at those bands and such models would be relevant. 

BuDAJ: No. However, this would be indeed very interesting. Thank you for the sugges- 
tion /information. 

RuciNSKi: I congratulate you for making order with definitions of albedo, as applied to 
stars and planets. The old definition was directly related to the so-called reflection effect 
in eclipsing binaries. 

Budaj : Thank you. 

SouTHWORTH: You find that the difference in radius between a spherical planet and the 

substcUar point of a Rochc-modcl planet is as large as 15%. When a planet is transiting 
we do not see its substellar point but we do measure its equatorial/ polar radius. This 
will be closer to the radius of a spherical planet so the problem is not as bad. Can you 
put a figure on this? 

BuDAJ: Yes. That is true (if one deals with transits and not with other phases). 
Rside/ Rpoie is not as big. I do not remember the exact numbers but they are in my 
paper for all the transiting planets. This ratio might be about 3% in the worst case. 
The effective radius of the planet (radius of the sphere with the same volume) is also 
tabulated there. 



